Anisotropic admixture in color-superconducting quark matter 
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The analysis of color-superconducting two-flavor deconfined quark matter at moderate densities 
is extended to include a particular spin-1 Cooper pairing of those quarks which do not partici- 
pate in the standard spin-0 diquark condensate, (i) The relativistic spin-1 gap A' implies sponta- 
neous breakdown of rotation invariance manifested in the form of the quasi-fermion dispersion law. 
(ii) The critical temperature of the anisotropic component is approximately given by the relation 
T' c ~ A'(T = 0)/3. (iii) For massless fermions the gas of anisotropic Bogolyubov-Valatin quasi- 
quarks becomes effectively gapless and two-dimensional. Consequently, its specific heat depends 
quadratically on temperature, (iv) All collective Nambu-Goldstone excitations of the anisotropic 
phase have a linear dispersion law and the whole system remains a superfluid. (v) The system 
exhibits an electromagnetic Meissner effect. 

PACS numbers: 12.39.Ki, 12.38.Aw,11.30.Qc 



Recent investigations suggest that the phase structure 
of QCD is very rich 0, Q • At low temperatures and high 
densities strongly interacting matter is expected to be a 
color superconductor At asymptotically high densi- 
ties, where the QCD coupling constant becomes small, 
this can be analyzed starting from first principles 0, , 
whereas at more moderate densities, present (presum- 
ably) in the interiors of neutron stars, these methods 
are no longer justified. In this region the low-energy 
dynamics of deconfined quark matter is often studied 
employing effective Lagrangians C e ff which contain local 
or non-local four-fermion interactions, most importantly 
interactions derived from instantons or on a more phe- 
nomenological basis 0, 0, 0- The non-confining gluon 
SU(3) C gauge fields are then treated as weak external 
perturbations, and neglected in lowest approximation. 

In this letter we consider the case of two flavors which 
is most likely relevant at chemical potentials just above 
the deconfinement phase transition. On physical grounds 
it is then natural to assume that C e g favors the spon- 
taneous formation of spin-0 isospin singlet Cooper pair 
condensates Q, 6 = (V> T C 75 t 2 A2 ip) , where ip is 
a quark field, C the matrix of charge conjugation, r 2 a 
Pauli matrix which acts in flavor space, and A2 a Gell- 
Mann matrix which acts in color space. Due to the latter 
SU(3) C is broken down to SU (2) c . This has the following 
consequences for the physical excitations of the system: 

(i) Corresponding to the mixing of the colors 1 and 
2 there are two Bogolyubov-Valatin quasiquarks for 
each flavor with the dispersion law Ef(p) = Ei^(p) = 
E ± (p) = yj(e p ±n) 2 + \A\ 2 . The energy gap A is the 
solution of a selfconsistent gap equation and is found to 
be typically of the order ~ 100 MeV in model calculations 
0, 0, |i| . e p = a/p 2 + M 2 , where M is an effective Dirac 
mass, related to the chiral condensate (ipip) via a selfcon- 



sistency equation |8( . For each flavor there is an unpaired 
quark of color 3 with the dispersion law e^(p) — e p ± 

(ii) Because of the spontaneous breaking of SU(3) C 
down to SU (2) c five of the eight gluons receive a mass 
(Meissner effect), whereas three remain massless [loj . 
Since no global symmetry is spontaneously broken there 
are no massless Goldstone bosons. 

(iii) The condensate 5 is invariant under a local U(l) 
transformation generated by Q = Q — t^tjAs, where Q is 
the electromagnetic charge operator and As a Gell-Mann 
matrix in color space. As long as this symmetry is not 
broken by other condensates, there is a "new" photon 
(a linear combination of the "normal" photon and the 
eighth gluon) which remains massless. This means, there 
is no electromagnetic Meissner effect. 

According to Cooper's theorem any attractive interac- 
tion leads to an instability at the Fermi surface. It is 
therefore rather unlikely, that the Fermi sea of color-3 
quarks stays intact. As only quarks of a single color are 
involved, the pairing must take place in a channel which 
is symmetric in color. Assuming s-wave condensation in 
an isospin-singlct channel, a possible candidate is a spin- 
1 condensate |y] . Although the size of the corresponding 
gap was estimated to be much smaller than A |fj , its ex- 
istence can have important astrophysical consequences. 
For example, if all quarks are gapped, the specific heat 
of a potential quark core of a neutron star (and hence 
the cooling of the star) is goverened by the size of the 
smallest gap Q. The same is true for other transport 
properties, like neutrino emissivity or viscosity. 

This letter is devoted to a quantitative analysis of this 
possibility. To this end we consider the condensate 

5' = (ip T Ca 03 r 2 P^iP) , (1) 
where a» v = i/2[Yn 1 '] and P 3 (c) = 1/3- 1/V3A 8 is 
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the projector on color 3. 5' is a ground-state expectation 
value of a complex vector order parameter <jp n = <f> n de- 
scribing spin 1 and breaking spontaneously the rotational 
invariance of the system. There are well-known examples 
for spin-1 pairing in condensed matter physics, e.g., su- 
perfluid 3 He, where some phases are also anisotropic [llj . 
In relativistic systems this is certainly not a very frequent 
phenomenon. It is possible only at finite chemical po- 
tential, which itself breaks Lorentz invariance explicitly. 
(Relativistic Cooper pairing into spin-1 with nonzero an- 
gular momentum was considered elsewhere, e.g., 0,^3-) 
Another example of non-isotropic quark pairing are crys- 
talline color superconductors [l3j. The role of spin-1 con- 
densates as an alternative to crystals in single color or 
single flavor systems has also been discussed in Ref . ■ 

Since rotational invariance is a global symmetry of 
the primary Lagrangian, an arbitrary small gap of the 
anisotropic phase implies specific gapless collective exci- 
tations with given Landau critical velocity crucial for the 
superfluid behavior of the system. We will briefly discuss 
this at the end of this Letter. It is also interesting to note 
that 6' is not neutral with respect to the "rotated" elec- 
tric charge Q and there is no generalized electric charge 
for which both, S and 6' , are neutral. This means, if both, 
5 and 5', were present in a neutron star, there would be an 
electromagnetic Meissner effect, which would strongly in- 
fluence the magnetic field. Recently, similar effects have 
been discussed in Ref. The detailed evaluation of 

the Meissner masses for our case remains to be done. 

For the quantitative analysis we have to specify the in- 
teraction. Guided by the structure of instanton-induced 
interactions (see, e.g , 0) we consider a quark-antiquark 
term 

C qq = G{(nf - (V^) 2 - (^ 7 5^) 2 + (^7 5 tV0 2 } (2) 
and a quark-quark term 

C qq = ~H S ^OCt 2 \ a 4, t )^ t COt 2 \ a ^ t ) 

0=75,1 

- HSa^CT 2 Xs^ T )(^ T Ca^T 2 \s^ T ) , (3) 

where and As are the antisymmetric and symmet- 
ric color generators, respectively. For instanton induced 
interactions the coupling constants fulfill the relation 
G : H s : H t = 1 : j : jg, but for the moment we will 
treat them as arbitrary parameters. As long as they stay 
positive, the interaction is attractive in the channels giv- 
ing rise to the diquark condensates 5 and 6' as well as to 
the chiral condensate (ify) . It is straight forward to cal- 
culate the mean- field thermodynamic potential Q(T,/i) 
in the presence of these condensates: 
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FIG. 1: M (dotted), A (dashed), and A' (solid) at T = as 
functions of the quark chemical potential fi using parameter 
set 1 (see text). The dashed-dotted line indicates the result 
for A' taking parameter set 2. 

where m is the bare quark mass, M = m — 2G{ijnl)), 
A = -2H S 5, and A' = AH t S'. The dispersion law for 
quarks of color 3 reads 

Ef(p) = \j^Ml ff +^T^ffY + \Kff\ 2 > (5) 

where fi 2 ^ — fi 2 + |A'| 2 sin 2 #, M e g = M/i/n e ff, and 
\A' eff \ 2 = |A'| 2 (cos 2 + M 2 /n 2 eff sin 2 6) . Here cosfl = 
P3/\p\. Thus, as expected, for A' ^ 0, E 3 (p) is an 
anisotropic function of p, clearly exhibiting the sponta- 
neous breakdown of rotational invariance. For M = 0, 
the gap A' e g vanishes at 9 = tt/2. In general its minimal 

value is given by A' = M\A'\/^J^ 2 + |A'| 2 . Expand- 
ing around its minimum the low-lying quasiparticle 
spectrum takes the form 

E3 (p±_ , pa) « ^ A' 2 + v 2 ± (p ± - Po ) 2 + v 2 p 2 , (6) 

where v x = (1 - (juA' 2 / {MA' 2 )) 2 ) 1 / 2 , v 3 = Po = 
^|A'|, and p\ = p\ + p\. This leads to a density of 
states linear in energy: 

1 u 2 + IA'1 2 
N(E) = - f 1 E 0(E - A' ) . (7) 

The actual values for A, A' and M follow from the condi- 
tion that the stable solutions correspond to the absolute 
minimum of f2 with respect to these quantities. Imposing 
dVl/dA'* — leads to the following gap equation for A': 

A' = 16tf t A'V /-AA^ (l±El)J T tanh^ , (8) 
J (2tt) 3 V s ' Ef 2T ' W 

where s = fi e ff(p 2 + M 2 g) x l 2 . Similarly one can de- 
rive gap equations for A and M by the requirements 
dQ,/dA* = and dfl/dM = 0, respectively. Together 
with Eq. |JH}, they form a set of three coupled equations, 
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FIG. 2: Ai/Ai(T=0) as function of T/Ai(T=0). Dashed: 
A, = A. Solid: Ai — A'. The calculations have been per- 
formed at fi = 450 MeV for parameter set 1. 



FIG. 3: Specific heat for parameter set 2 at /j, = 450 MeV 
as function of T/T^.. Solid: full calculation, dashed: result 
for M = 0, dotted: without spin-1 condensate. The dashed- 
dotted line indicates the result of Eq. (1101 . 



which have to be solved simultaneously. However, the 
equations for A and A' are not directly coupled, but 
only through their dependence on M. 

In our numerical calculations we use a sharp 3- 
momentum cutoff A to regularize the integrals. We then 
have five parameters: m, A, G, H s , and H t . We choose 
m = 5 MeV, A = 600 MeV, and GA 2 = 2.4 - lead- 
ing to reasonable vacuum properties, M = 393 MeV and 
(uu) = (— 244MeV) 3 -, and the instanton relation to fix 
H s and H t ("parameter set 1"). The resulting values of 
M, A, and A' as functions of /i at T = are displayed 
in Fig. n The chemical potentials correspond to baryon 
densities of about 4-7 times nuclear matter density. In 
agreement with earlier expectations |6j A' is small com- 
pared with A. However, its value rises strongly with /i. 
Being a solution of a selfconsistency problem, A' is also 
extremely sensitive to the coupling constant H t . If we 
double the value of H t ( "parameter set 2" ) , we arrive at 
the dashed-dotted line for A', which is then comparable 
to A. As a consequence of the factor (1 — Pj_/s) in the 
gap equation |JSJ, A' is very sensitive to value and the 
form of the cutoff. 

With increasing temperature both condensates, 6 and 
5', are reduced and eventually vanish in second-order 
phase transitions at critical temperatures T c and T' c , re- 
spectively. It has been shown [j| that T c is approximately 
given by the well-known BCS relation T c ~ 0.57A(T=0). 
In order to derive a similar relation for T' c we inspect the 
gap equation © at T = and in the limit T —> T' c . Ne- 
glecting M (since M <C /U this is valid up to higher orders 
in M 2 //it 2 ) and antiparticle contributions one gets 
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Since the integrand is strongly peaked near the Fermi 
surface, the |p| -integrand must approximately vanish at 
\p\ = fi, after the angular integration has been performed. 



From this condition one finds to lowest order in A'///: 
T'JA'{T = Q) « 1/3 . The analogous steps would lead 
to T c /A(T=0) w 1/2 instead of the textbook value of 
0.57. This gives a rough idea about the quality of the 
approximation. Note that there are other ex amp les of 
diquark condensates, where T c ^ 0.57A(T=0) 16|. This 
is also the case for crystalline superconductors [17|. 

Numerical results for A(T) and A'(T) are shown in 
Fig. [2 The quantities have been rescaled in order to 
facilitate a comparison with the above relations for T c 
and T' c . Our results are in reasonable agreement with 
our estimates. These findings turn out to be insensitive 
to the actual choice of parameters. 

The specific heat is given by c v = -Td 2 Vt/dT 2 0. 
For T -C T c it is completely dominated by quarks of color 
3, since the contribution of the first two colors is sup- 
pressed by a factor e~ A / T . Neglecting the T-dependence 
of M and A', and employing Eq. J7|), one finds 
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which should be valid for T <C T' c . In this regime c„ 
depends quadratically on T for T > Ag, and is expo- 
nentially suppressed at lower temperatures. To test this 
relation we evaluate c v (T) explicitly using Eq. JIJ . The 
results for fixed /i =450 MeV are displayed in Fig. EI F° r 
numerical convenience we choose parameter set 2, leading 
to a relatively large A'(T=0) = 30.8 MeV. The critical 
temperature is T' c ~ 0.40 A'(T=0). For the energy gap 
we find A' = 0.074 T' c . It turns out that Eq. evalu- 
ated with constant values of A' and M, (dashed-dotted) 
is in almost perfect agreement with the numerical result 
(solid) up to T m T' c /2. The phase transition, causing 
the discontinuity of c„ at T = T' c , is of course outside 
the range of validity of Eq. (|10|) . We also display c v for 
M = (dashed). Since A vanishes in this case there is 
no exponential suppression, and c v is proportional to T 2 
down to arbitrarily low T . However, even when M is in- 
cluded, the exponential suppression is partially cancelled 
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by the sum on the r.h.s. of Eq. I|10|) . For comparison we 
also show c v for a system with A' = 0, which exhibits a 
linear T dependence at low temperatures (dotted). 

Our results show that, even though the magnitude of 
the gap parameter A' is strongly model dependent its re- 
lations to the critical temperature and the specific heat 
are quite robust. Thus, if we had empirical data, e.g., 
for the specific heat of dense quark matter, they could 
be used to extract information about the existence and 
the size of A'. In this context neutron stars and their 
cooling properties are the natural candidates to look at. 
In Ref. M it was suggested that the exponential suppres- 
sion of c v related to the potential pairing of quarks of 
color 3 might have observable consequences for the neu- 
trino emission of a neutron star. This argument has to 
be somewhat refined since, as seen above, c v (T) first be- 
haves as T 2 and the exponential suppression sets in only 
at T < Aq. The relevance of c v and the possible ef- 
fect of diquark condensates on neutron star cooling was 
also discussed in Ref. Ilfll . On the other hand it has 
recently been argued |2fl |. that the constraints imposed 
by charge and color neutrality might completely prohibit 
the existence of two-flavor color-superconducting matter 
in neutron stars. 

Because of the spontaneously broken U(l) x 0(3) sym- 
metry in Eq. (JJJ, for A' ^ there should be collective 
Nambu-Goldstone excitations in the spectrum. However, 
due to the Lorentz non-invariance of the system there can 
be subtleties 0, H3, . The NG spectrum can be ana- 
lyzed within an underlying effective Higgs potential 

V(</)) = -a 2 4^n + ^A 1 (4^ n ) 2 + -A 2 44^ m ^ ro , (11) 

for the complex order parameter <f> n , with \i + X 2 > 
for stability. For A2 < the ground state is character- 
ized by (fivac = (f^) 1 ^ 2 (0, 0, 1) which corresponds to our 
ansatz Eq. for the BCS-type diquark condensate 5'. 
This solution has the property (S) 2 = (<$alS(j>v}c) 2 = 0. 
The spectrum of small oscillations above (j^vac consists 
of 1+2 NG bosons, all with linear dispersion law: one 
zero-sound phonon and two spin waves |22^ . Implying a 
finite Landau critical velocity, this fact is crucial for a 
macroscopic superfluid behavior of the system [23| . 
Note that for A2 > there is a different solution 
L, x a l\ 1/2 (M,0) with (S) 2 = 1. In this case 



6 ( 2 ) 
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the NG spectrum above <jr v dc consists of one phonon with 
linear dispersion law and one spin wave whose energy 

tends to zero with momentum squared |22^ . The quasi- 

(2) 

quark dispersion law corresponding to cpvdc has recently 
been discussed in Ref. 14]. A detailed analysis of the 
transport properties implied by this type of pairing seems 
of particular interest for the phenomenology of neutron 
stars since the dispersion law allows for gapless excita- 
tions even for massive quarks (l^j . 
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